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While there are many forms of time-transfer, the most precise long-distance forms currently used for the 
generation of International Atomic Time (TAI) and Coordinated Universal Time (UTC) involve either GPS 
or Two Way Satellite Time and Frequency Transfer (TWSTFT). This paper gives a brief description of 
their current and future capabilities, with emphasis on how their uncertainties affect UTC. Some of these 
uncertainties are due to inherent modeling or receiver instabilities, while others can be reduced through 
temperature and humidity stabilization, electronic impedance matching, and multipath minimization. The 
residual time-transfer uncertainties directly affect the uncertainties in each individual laboratory’s 
realization of UTC. 

GPS is the best-known means of precise time transfer, and the introduction of Common-View GPS (CV) 
contributed significantly to the great increase in the stability of TAI over the last 15 years. GPS 
pseudorange observations are currently adjusted by correcting the broadcast orbital and ionosphere delay 
information with post-processed values, and reduced in common-view mode. However “melting-pot” 
techniques have been shown to be more precise over long baselines and “P3” observations, using dual¬ 
frequency carrier-phase receivers, allow use of measured ionosphere delays. P3-links have been used by 
the International Bureau of Weights and Meaures (BIPM) for TAI-generation; plans to greatly increase 
their role have been circulated. The continuous GPS observations are supplemented by periodic receiver 
calibrations, which have accuracies of several ns. Greater precision (20 ps), but not significantly greater 
accuracy, would be attainable if carrier-phase GPS techniques were employed [1]. Several institutions, and 
in particular the International GNSS Service (IGS), routinely provide publicly available carrier-phase based 
timing solutions in near-real time [2], Other institutions will reduce any laboratory’s data for free via an 
automated procedure. A few institutions, such as JPL/NASA and NRCan, provide real-time carrier-phase 
based time transfer as well [3,4] 

Beginning in 2000, time-transfer links using TWSTFT replaced some GPS links as the primary operational 
link, and currently over half the clocks used for TAI-generation are linked to other sites via a direct 
TWSTFT link. TWSTFT links can be calibrated to achieve slightly subnanosecond absolute time transfer, 
with a precision in the range of a few hundred ps [5]. Carrier-phase TWSTFT holds the promise of ps-level 
precision, but this promise has not yet been realized in any system [6]. 

The generation of TAI involves merging information via dissimilar links, which complicates the relevant 
statistical approach and increases the dependence of TAI upon the stability of the time-transfer equipment 
at those “pivot sites”, which are linked to some laboratories via TWSTFT but linked to other laboratories 
via GPS. The data from the links are used to transfer the time to generate UTC, which expressed as UTC - 
UTC (k). The first step in this process is to generate the free atomic time scale (EAL). EAL is defined 
using the ALGOS algorithm [7-9] as: 

EAL{t ) = ^ w’ \h t (t) + /z/(r)] (1) 


where N is the number of the atomic clocks, w, the weight of the clocks, h^t) is the reading of clock H, at 
time t, and V (t) is the prediction of the reading of clock //, to guarantee the continuity of the time scale. 
The weight attributed to a given clock reflects its long-term stability, since the objective is to obtain a 
weighted average that is more stable in the long term than any of the contributing elements [10]. The 
ALGOS software package is used in the Time Section of the BIPM in order to generate EAL, and 
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afterwards UTC is generated by the addition of leap seconds and frequency steers. The data used by 
ALGOS take the form of the time differences between readings of clocks, written as: 


x iJ (t)=h j (t)-h i (t). (2) 

This results in a system of N equations and N unknowns, the solution yielding: 

Xj (/) = EAL(t) -hj=Yj Wf \h[ (t) - X'J (f )] (3) 


If we choose a particular clock H,, we can see that the difference between that clock and EAL depends on 
weights, clock prediction, and measured clock differences. The predictions and the weights are fixed by 
appropriate algorithms based on the past clock behavior and therefore can be considered as time-varying 
deterministic parameters. Suboptimal estimation of these parameters would affect the uncertainty of TAI as 
realization of the Terrestrial Time ( TT ), but they do not affect the knowledge of the difference between 
EAL and clock Hj. The measures x kJ are thus the only contributors to the uncertainties in Xj. We ignore the 
contribution of the uncertainty given by measurements of clocks located inside the same laboratory, since 
they are both negligible and indistinguishable from the noise of the individual clocks. The generation of 
UTC from EAL through the addition of pre-determined leap seconds and frequency steers does not add to 
the uncertainty. The uncertainties of [UTC - UTC (k)] are therefore identical to the uncertainties of [TAI- 
UTC (k)] and [EAL - UTC (k)], and we conclude that the uncertainties of the links among laboratories are 
the only source of the uncertainty of UTC - UTC (k). 

In order to correctly model the processes, the correlation of time-transfer links that share common sites 
must be estimated and accounted for. One approach is to use the law of propagation of uncertainties 
[11,12]. The uncertainty in the xf t) can be found for a generic quantity y measured by means of direct 
measurements of the input quantity x„: 

y = f(x x,x 2 ,....,x m ). (4) 


The expression of the law of the propagation of uncertainty [12] is given by: 

8f 8f 

«■<«> (5) 




where the first term corresponds to the effect of the uncertainties of the input quantities x„ and the second 
term accounts for the correlation between them. 


In our case, we are interested in evaluating the uncertainty of the quantity x k = [EAL - UTC (k)] which here 
plays the role of the indirect quantity y, and the uncertainty contributions are only due to the measurement 
noise of the links x hJ (t): 




► f -±*K + 


where u], = u 2 EAL _ h , . 

The weights of the clocks are available from the BIPM Web site, and the uncertainty of links between the 
clocks [13] are published in Circular T (see the BIPM Time Section’s FTP server at www.bipm.org ). It can 
be demonstrated that the ALGOS algorithm would generate the same results if each laboratory’s clocks 
were replaced by a single “equivalent” clock whose reading was the weighted average of the individual 



clocks and whose weight in EAL was the sum of the individual clock weights. However, even if the 
uncertainties of transferring time internally within the labs were not considered part of the noise of the 
individual clocks, the formula could still be applied, as the double summation would account for the 100% 
correlation of the time transfer noise between clocks in the same lab. 

For application to the computation of UTC, we need to consider the correlations between the links. If all 
time transfer were achieved using a single system per site, and if all sources of noise were site-based, such 
as mostly happens for Melting Pot GPS, also termed All-in-View (AV) [14,15], then all possible links 
would obey the following closure relation: 

x ij (t) + x jk (t) + x ki (t) = 0 (7) 

In this situation, the noise of each site’s time transfer system would be indistinguishable from the noise of 
its clocks. The dominant uncertainty would be given by the site noise and all the external clocks would be 
seen through that dominant noise independently of their location. This noise would affect UTC -UTC (k) 
itself according to the total weight of the site clocks, but due to the closure relation it would not affect the 
difference between the clocks of any other two sites. Therefore, in such a situation: 

U UTC-UTC(k) = (1 — W k ) U k,k (8) 

where is the uncertainty of site k’s time transfer system. 

Uncertainties in time transfer using GPS common view (CV) are largely, though not entirely, laboratory- 
based. However, even for CV observations made using every available satellite, closure violations will 
arise if simultaneous satellite observations are recorded at only two of the three sites. In such cases, orbit 
mis-estimation and receiver noise will contribute uncertainties, and any azimuth or elevation-dependent 
asymmetries in the multipath environment would cause both uncertainties and biases. Since calibration is 
achieved by an all-sky sampling that is systematically different from the sky-sampling of CV, systematic 
multipath will also lead to uncertainties. Despite these noise sources, the closure relation largely holds for 
common view, and the largest source of uncertainty is typically due to variations of the receiver system that 
are common to all the data. In a BIPM-sponsored study [15] of 40 triplets of European and North 
American sites, the median closure violation was about 1 ns RMS for CV, and only 350 ps RMS for AV. 

For Two-Way Satellite Time and Frequency Transfer (TWSTFT), the noise is again largely site-dependent. 
Some closure violations can occur because the observations between pairs of sites are typically made at 
different spread-spectrum codes and slightly different frequencies. The largest source of closure errors is 
probably due to the fact that the received signals are shaped by the product of the transmitting and 
receiving bandpasses, while the delay and certain noise components such as the cable-dependent multipath 
can systematically vary over the bandpass [16]. However, such closure violations can be reduced somewhat 
through baseline-dependent calibrations. In those cases where a TWSTFT system is calibrated with GPS, 
the uncertainties in the calibration are determined by the uncertainties in the GPS calibration. For 
calibrated data, the median closure violation of 434 European and N. American site-triplets is 680 ps RMS. 

Special situations arise when one site is a pivot site, connected to some sites by one technique and other 
sites by a different technique. Let us assume a bias B exists in the GPS equipment at the pivot site. In this 
case, it is easy to show that the bias would affect those laboratories linked to the pivot site through the GPS 
system as follows: 

A UTC - UTC(k) = (1 — W G )B (9) 

where W G is the sum of all the weights of the laboratories linked to the pivot site by its GPS system. 

The bias would affect the pivot laboratory and those linked to it by the pivot’s TWSTFT system as follows: 

A UTC -UTC(k) = —W G B (10) 



Under the normal circumstances described in this paper, the existence of any biases would not carry any 
significant statistical implications, as they would be directly related to the tabulated uncertainties in the 
links themselves. However, the above equations illustrate the dependence of TAI and UTC upon the 
equipment at any such “crossover pivot”, and in the case of equipment failure or aging, such dependencies 
can result in non-Gaussian behavior. 

In order to optimize the use of time-transfer in TAI-generation, several algorithms are under consideration. 
One algorithm would be to use all available links, rather than the minimum required to generate TAI, in a 
large least squares fit in which bias-parameters would be associated with all crossover sites. The bias 
parameters represent the bias between the site’s different systems (typically GPS and TWSTFT, but they 
could also be the difference between the SPS GPS receiver and a geodetic receiver), and the weighted 
average of those bias parameters would be constrained to be zero. The calibration of all links would be 
pre-set by a self-consistent methodology. More elaborate schemes are also possible, but the common 
theme is that TAI-generation can be improved through improved modelling and removal of the systematic 
errors. As our frequency standards become more precise and more accurate, such techniques will become 
more necessary. 
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